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FOREWORD 


This is Volume II of three separately bound volumes in which are reported the 
research completed under the general terms of the Office of Civil Defense Subtask 
Number 4113E, "Sensitivity Analysis of Civil Defense Systems and Components." 

The author acknowledges the valuable assistance of Mr. Philip McMullan in 
completing the research and of Mr. McMullan and Mr. Robert Brooks in preparation 
of the final report. 






ABSTRACT 


In order to perform a sensitivity analysis of selected parameters of interest 
in Civil Defense systems analysis, probable casualties are estimated for 8 SMSA's 
over a range of fallout environments and shelter utilization patterns. The selected 
parameters are: SMSA population, population density, and ratio of shelter spaces to 
population; fallout arrival time and reference intensity; and restrictions on move¬ 
ment of people to shelter, leading to varying patterns of shelter utilization. The 
SMSA's are selected by "judgment sampling" and range in population from 74,000 to 
408,000. The fallout environments used range from a reference intensity of 600 r/hr 
and 7 hours time of arrival to a reference intensity of 30,000 r/hr and 1 hour time 
of arrival. The movement-to-shelter restrictions are: (1) movement restricted to 
the Standard Location (SL) of residence, (2) movement restricted to within two miles 
of the SL of residence, and (3) unrestricted movement to shelter anywhere within the 
SMSA. Also, (4) the transportation algorithm is used to determine the optimal 
(minimum casualty) allocation of people to shelter for each time of arrival and 
reference intensity combination. This allocation serves as a benchmark of ideality 
against which to measure other patterns of shelter utilization. Casualties are 
computed for each of the four movement patterns over the range of attack environ¬ 
ments. It is concluded that efficient shelter utilization is very important when 
shelter spaces exceed population; inefficient utilization can lead to a large number 
of avoidable casualties. Movement through fallout for a short time in order to 
reach a better shelter theoretically is warranted in many circumstances. 
































LIST OF TABLES 


Table 

I Generalized Matrix for ERD "Costs" 


Page 


3 


II Residential Populati on and Fallout Shelter Characteristics 

for Eight SMSA's . 6 

III Fallout Environments Used in 8 SMSA Analyses . .... 10 

IV Expected Casual tier for Each Fallout Environment in 

Orlando, Florida ..... 13 

V Expected Casualties for Each Fallout Environment in 

Tucson, Arizona . ... IB 

VI Expected Casualties for Each Fallout Environment in 

Charleston. South Carolina . ... 22 

VII Expected Casualties for Each Fallout Environment in 

Mobile. Alabama . 26 

VIII Expected Casualties for Each Fallout Environment in 

Macon, Georgia . ..... 30 

IX Expected Casualties for Each Fallout Environment in 

Pittsfield. Massachusetts . 34 

X Expected Casualties for Each Fallout Environment in 

Waterloo. Iowa . ....... 36 

XI Expected Casualties for Each Fallout Environment in 

Richmond. Virginia . 41 

XII Summary of Casualty Estimates ( 7 . of Population) for Reference 

Intensity of 2000 r/hr and Time of Arrival in 1.5 Hours ... 48 




















LIST OF FIGURES 


Figure Page 

1 Distribution of Population Density .. . 5 

2 Distribution of Shelter/Population Ratio . . 5 

3 Joint Distribution of Population Density and Shelter/Ponulation 

Ratio of 200 SMSA’s. 7 


4 _Cumulative Distribution Function for Reference Intensity . . .. 8 

5 Cumulative Distribution Function for Time of Arrival of 

Fallout . 9 


6 Cumulative Distribution of Population and Shelter in 

Orlando, Florida . . . ....... 12 

7 Percentage of Casualties for Each Fallout Environment in 

Orlando, Florida. 14 

8 Percentage of Population for Each PF Category in Orlando, 

Florida. 14 

9 Cumulative Distribution of Population and Shelter in 

Tucson, Arizona. 15 

10 Percentage of Population Sheltered for Each PF Category 

in Tucson, Arizona . ................. 16 

11 Percentage of Casualties for Each Fal out Environment in 

Tucson, Arizona .. ^. 17 

12 Cumulative Distribution of Population land Shelter in 

Charleston, South Carolina . . . . . . . 19 


13 Percentage of Population Sheltered for Each PF Category in 

Charleston, South Carolina ........... .. . 20 

14 Percentage of Casualties for Each Fallout Environment in 

Charleston, South Carolina.... 21 

15 Cumulative Distribution of Population and Shelter in Mobile, 

Alabama.... 23 

16 . Percentage of Population Sheltered for Each PF Category in 

Mobile, Alabama . . ........ 24 


xi 



















LIST OF FIGURES (CONTINUED) 


Figure 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 
B-l 


Page 


Percentage of Casualties for Each Fallout Environment in 

Mobile, Alabama ...... . . . 25 

Cumulative Distribution of Population and Shelter in 

Macon, Georgia . ..... 27 

Percentage of Population Sheltered for Each PF Category in 

Macon, Georgia ..... . . 28 

Percentage of Casualties for Each Fallout Environment in .. 

Macon, Georgia . 29 

Cumulative Distribution of Population and Shelter in 

Pittsfield, Massachusetts . . . 31 

Percentage of Population Sheltered for Each PF Category in 
Pittsfield, Massachusetts ....... . ... 32 

Percentage of Casualties for Each Fallout Environment in 
Pittsfield, Massachusetts . 33 

Cumulative Distribution of Population and Shelter in 

Waterloo, Iowa. 35 

Sensitivity of Casualty Calculations to Movement Resirictions-- 
Waterloo, Iowa. 37 

Percentage of Population Sheltered for Each PF Category in 
Waterloo, Iowa . .. 38 

Percentage of Casualties for Each Fallout Environment in 
Waterloo, Iowa .. 39 

Cumulative Distribution of Population and Shelter in Richmond, 
Virginia. 40 

Percentage of Population Sheltered for Each PF Category in 
Richmond, Virginia. 42 

Percentage of Casualties for Each Fallout Environment in 
Richmond, Virginia ........ . ........ 43 

Trade-Off Between PF and Time in Fallout.. 45 

Trade-Off Between PF and Fallout Environment . . 46 


Illustration of Finite Difference Formulation . . . B-l 























A Sensitivity Analysis of Selected Parameters 
Based on 8 SMSA's 


I. INTRODUCTION 

The purpose of this volume is to analyze the sensitivity of estimates of the 
effectiveness of fallout shelter systems to changes in movement assumption? and 
resultant shelter utilization under various fallout environments, shelter/population 
ratios, and population densities. 

Eight Standard Metropolitan Statistical Areas (SMSA's) were chosen representing 
a range of population densities and shelter/population ratios. Four modes of alloca¬ 
ting population to shelters, based or, travel restrictions in time or distances, and 
six types of fallout environment were assumed for each of the SMSA's. Expected 
casualties were computed for each combination of allocation mode and fallout environ¬ 
ment. 

The results were then analyzed to answer these and other questions: 

(1) How effective are various pat*?rns of shelter utilization in each of the 
SMSA's? 

(2) What trade-offs exist between protection factors and movement in fallout? 

II. DEFINING THE MODEL 

The four types of movement restrictions considered in the analysis are as 
follows: 

(1) Movement restricted to within each standard location (SL) and no movement 
in fallout necessary to reach shelter. 

(2) Movement restricted to within an SL or to SL's not more than two miles 
(distances are rectangular distances--rather than straight line--measured 
between the centers of SL's); no restriction on speed of movement. 

(3) Movement unrestricted within the SMSA and no restriction on speed of move¬ 
ment. 
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(4) Movement unrestricted within the SMSA and a constant movement speed of two 
miles per hour beginning with weapon detonation; movement in fallout per* 
mitted; casualties minimized using the transportation algorithm (see below). 

Assumptions for all cases were as follows: 

(1) All standard locations in an SMSA have the same fallout environment; 

(2) The fallout intensity increases linearly during the time fallout particles 
are being deposited; 

(3) Maximum Equivalent Residual Dose (EPD) is calculated by a finite differencing 
method which examines dose accumulation and recovery from time of arrival 
through the initial shelter period (see Appendix B); 

(4) Midpoint values of each Protection Factor (PF) category are used in casualty 
calculations; 

(5) A person experiencing a Maximum ERD of 200r is considered a casualty. 

Calculations of casualties for the case with no movement outside the standard 
location were straightforward. NFSS Phase 1 data [Reference 1] were used for each 
SL to determine the population and shelter spaces of each PF category, including 
Category 1. The population was then allocated manually on the basis of utilizing 
the highest available PF shelter space. Any population in excess of the number of 
shelter spaces available was considered to have a PF of 2. 

In the third type of movement restriction (movement anywhere within the SMSA 
and no speed restriction), population was allocated to the highest available PF 
shelter spaces, regardless of their initial locations within an SMSA. 

In the case of movement restrictions of two miles, it was necessary to compute 
the distances between centers of SL's. The rectangular distance was used since it 
is reasonable that the route to shelters in most cities must follow some right angle 
street pattern rather than a straight line. The allocation of population to shelters 
was then done manually, always using all the space with the highest PF within the 
two-mile movement restriction. 

The fourth case of movement restriction, limiting movement speed to two miles 
per hour, was much more complex. It was necessary to include the possibility of 
movement in fallout in order to reach a more effective shelter. In order to cal¬ 
culate the movement pattern leading to the minimum number of casualties, it was 
first necessary to calculate the ERD accumulated in moving from one standard loca¬ 
tion (S^) to shelter in another standard location (Sj) of protection factor (It). 











When the travel tine between standard locations was great enough to require movement 
in fallout, it was necessary to determine this exposure time and to calculate the 
dose sustained outside, as well as that sustained while in shelter. 

The process of determining the EPD for movement between standard locations under 
various fallout environments was programed for the CDC 3600 computer. The generalized 
matrix that was used for these ERD "costs'* is shown below in Table I. 


TABLE I 

Generalized Matrix for ERD "Costs" 



S 1 


S 2 


Destinations 

s j 


Origins 

PF 1 

pf 2 

PF l 

pf 2 

• • • • •» PF^ 

• • • • 

PF k 

S 1 

S lll 

S 112 

S 121 

S 122 

.. .. .. Sj 

• • • • 

s i Jk 

S 2 

• • 

S 211 

• • 

S 212 

• • 

S 221 

S 222 

. S 2jl 

• • • • 

S 2jk 

• • 

• • 

• • 

(e.g., S ^22 
destination 2 

the"ERD" cost to move from origin 1 

and enter a PF category 2 space) 

to 

S i 

S ill 

S il2 

S i21 

S i22 

. S ijl 

• ♦ • • 

S. .. 
ijk 


With such a matrix, and with the assumption that casualties occur when 
ERD > 200r, the transportation problem method of linear programming was used to mini¬ 
mize the number of casualties within each SMSA for each fallout environment. The 
mathematical description of the transportation model is given in Appendix B. 


III. DATA FOR EVALUATION 



A. Selection of Local Areas (SMSA's) 


Eight SMSA's were selected according to two criteria, shelter/population ratio 
and population density, by a procedure known as judgment sampling [Reference 2], 
Data from judgment samples only suggest or indicate conclusions, which is precisely 
the purpose intended in this analysis. 

Prior to selecting the SMSA's, distribution functions of population density 
and shelter/population ratio were considered (Figures 1 and 2). The distributions 
represent 200 of the 213 SMSA's listed in the National Location Code [Reference 3]. 
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Fig. 1. Distribution of Population Density (200 SMSA's) 
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Shelter/Population Ratio (spaces/person) 


| Fig. 2. Distribution of Shelter/Population Ratio (200 SMSA's) 







Two constraints were placed on the sampling procedure: 

1. SMSA's with population exceeding 500,000 were excluded. This constraint 
was dictated by the time required to solve the algorithm used to allocate 
people to shelter. 

2. Multiple SMSA's (e.g., Greensboro-High Point, N. C.) were excluded. In 
most cases, these areas require special consideration because of complex 
movement restrictions. 

Hie shaded portion in Figures 1 and 2 shows the distribution (after elimination 
of multiple SMSA's and SMSA's with population greater than 500,000) frcm which the 
judgment sample was taken. The characteristics of the selected SMSA's are listed in 
Table II, Residential Population and Fallout Shelter Characteristics for Eight SMSA's . 

In order to visualize the characteristics of the selected SMSA's and their 
relationships, consider the joint distribution of population density and shelter/ 
population ratio as shown in Figure 3, Joint Distribution of Population Density and 
Shelter/Population Ratio for 200 SMSA's. If we consider the sample to represent all 
SMSA's lying within the boundaries formed by the line segments connecting each 
selected SMSA, then 146 of the 200 will be represented. (It should be pointed out, 
however, that all SMSA's within the boundaries described do not satisfy the popula¬ 
tion constraint of being less than 500,000; also some multiple SMSA’s are included.) 

B. Selection of Fallout Environments 

One form of available data providing reference intensity (I q ) and associated 
arrival times for fallout (T^) for all SMSA's is the RISK war game analysis 
[Reference 4], In order to use representative numbers for the fallout environment, 
data from an 0CD study in a form similar to that described in Reference [4] were used 
to derive distribution functions for both I and T. (Appendix A). Figures 4, Cumu,- 
lative Distribution Function for Reference Intensity, and 5, Cumulative Distribution 
Function for Time of Arrival of Fallout, relate the probability of occurrence at each_„. 
level of I q and T^, respectively, for the 0CD study data. The sample size (N) was 
20,000; i.e., 100 attacks on each of 200 SMSA's. The available data were aggregated 
in such a way that it would be inaccurate to associate a joint probability of occur¬ 
rence of I q and T^. However, we have assumed a lower limit on the probability of 
occurrence of a given combination of I q and T^ to be the product of their individual 
probabilities. By assuming complete dependence, we can arrive at some upper limit 
on the probability of occurrence for the same combination. Thus, we can describe 
some probability of ~ terval within which the joint probability of I Q and T^ must lie. 













People 
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Joint Distribution of Population Density and Shelter/Population Ratio for 200 SMSA's 
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Cumulative Distribution Function for Reference Intensity 












































































Time of Arrival, (hours) 

Fig. 5. Cumulative Distribution Function for Time of Arrival of Fallout 


















For example, the fallout environment will be I < 3000 r/hr and T, > 1 hour between 

O “ A “ 

81 percent and 90 percent of all occurrences (see Table III for method of computing). 

Six fallout environments (combinations of and I q ) were selected and used to 
test the effectiveness of various shelter systems (see Table III). 

TABLE III 

Fallout Environments Msed in 8 SMSA Analyses 


I T 

o A Approximate Range of 

(r/hr) (hr) Cumulative Probability 


< 600 


£5.5 


367* 

- 

** 

607. 

2000 

-.. • 

1.5 

-...-. 

727. 


857. 

3000 


1 


817. 

- 

907. 



*** 





5000 


1 


867. 

- 

957. 



*** 





10000 


1 


897. 


9S7. 

30000 


IHrk 

1 


.907. 

- 

99.97. 


Lower Bound ; Since it is generally true that late arrival times correspond to 
low values of I c »it approximately can be said that: 

P (t a > t | i o < i) > P (t a > t) 


so that 

P [I o < i, T A > t) - P (I o < i) • P {T a £ t | I Q < i) > P (I o < i) • P (T a > t) 
for example: 

P {I < 600, T. £ 5.5) > P (I < 600) • P {T. > 7} £ (.60) (.60) £ .36 

Upper Bound : Since a probability is always <1, it is true that: 

P [I o < i, T A > t) - P (I o < i) • P (T a > t [ I D < i) < P (I 0 <Tf *1 ; 
and similarly 

P (I o < i, t a > t) < 1 • P (T a > t) 
so that 

P {I 0 < 600, T a >5.5 }< min [p {I q < 600 ) , P (T A > 7 ] ] < .60 

In some cases P (T A > t) was not known. In this case, the uppe. bound 
P (I Q < i, T a > t) £ P (I q < i) should be used. 
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IV. SM$A ANALYSIS 


A. 


Casualty Computations and Shelter Utilization Assumptions 

The model of SMSA's described in the previous sections was used to calcul 
the number of casualties in each SMSA at each of the fallout environments desc 
in Table III. (The term "casualty" in this report refers to anyone receiving 
Maximum ERD greater than 200r.) Casualties were computed for each of the four 
ment and shelter allocation rchemes and were used to rank their relative effec 
ness. 


ate 

ribed 

a 

move- 

Itive- 


Only one of the four movement-to-shelter procedures allowed for movement 
fallout. In each of the four procedures, those persons" "notsheltered in NFSS 
(including PF Category 1) were gi^en a PF of 2. ^Residential basements were no 
employed in these analyses. 


Ln 

spaces 


It is recognized that residential basements afford adequate protection in 
intensity fallout environments. However, excluding these spaces will only affect 
the casualty computations in one fallout environment of our analysis; i.e., at 
reference intensity equal to 600r/hr and time of arrival equal to 5.5 hours, rtsiden 
tial basements afford adequate fallout protection. Three SMSA's of the sample 
shelter in residential basements for more than 107. of the population—Pittsfield, 
Waterloo, and Richmond [Reference 5]. Of the thre|, only Pittsfield has a shelter 
deficit. Therefore, in the analysis presented here, only in Pittsfield could the 
total number of casualties be reduced significantly by using residential basement 


low 


have 


data, and then only in the least severe fallout environment. 


B. 


Analysis of the 8 SMSA's 1 

The following paragraphs present a casualty analysis for each of the 8 SM£>A's. 


They are presented in ascending order of shelter spaces to population ratio. Casu¬ 
alties for the four movement assumptions and eight fallout environments are presented 
with and without spaces obtainable in NFSS shelter through added ventilation or 
shielding. An analysis is made for each SMSA of the effect on casualties of the 
movement assumptions. 


1. Orlando. Florida 


e four 


The SMSA has a residential population of 263,540 and a population density 

I 

of 4170 people per square mile. The population density is slightly less than 
the average of the 200 SMSA's. 


I 


- 11 











Fig. 6. Cumulative Distribution of Population and Shelter 
in Orlando, Florida 

Approximately 20 percent of the population can be sheltered in spaces 
located in Phase 1 of the NFSS. Only 5 percent of the population can be 
sheltered in spaces having a protection factor of at least 100 and minimum 
volume requirements (500 cubic feet per person with no ventilation or 10 square 
feet with adequate ventilation). __ 

The plot of the cumulative distribution of population and shelter as a 
function of radial distance from the city center shown in Figure 6 indicates 
a shelter deficit throughout the SMSA; further, almost no shelter spaces exist 
beyond a two-mile radius of the population center. Due to the small number of 
available shelter spaces, both the allocation restricting movement to two 
miles and the constant movement speed restriction result in the same number of 
casualties as that if movement were unrestricted. 

Casuilty computations are significantly affected by the assumption of 
movement restricted to within the SL; at least this is true for the less 
severe fallout environments. The expected number of casualties at each fallout 
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environment is shown in Table IV. Figure 8, Percentage of Casualties for Each 
Fallout Environment in Orlando, Florida, shows expected casualties from Table 
IV as a percentage of the total SMSA's population. 


TABLE IV 

Expected Casualties for Each Fallout 
Environment in Orlando. Florida 


Fallout 

Environment 

Expected Casualties (7.) 

Reference 

Arrival 

Movement 

Movement 

Constant 


Intensity 

Time 

Restricted 

Restricted — 

— Movement — • 

- - - Movement.. 

(r/hr) 

Oir) 

to own SL 

to 2 Mi. 

Speed 

Unrestricted 

600 

5.5 

94 

78 

78 

78 

2000 

1.5 

94 

78 

. 78 

78 

3000 

1 

96 

89 

89 

89 

5000 

1 

98 

93 

93 

93 

10000 

1 

98 

.98 

98 

98 

30000 

1 

99 

99 

99 

99 


Data from Phase II of the NFSS survey [Reference 6] points up that 
approximately 7 percent of the population could be sheltered ii improvable 
spaces. The cost of improvement would be $91,000, or approximately $6.50 per 
space. A graphical representation of the shelter status is presented in 
Figure 7, Percentage of Population Sheltered for Each PF Category in Orlando, 
Florida. 

___It is quite obvious that this area would suffer severely-in any fallout- 

environment. Also, only slight advantage would be gained by improving existing 
spaces. The only solution appears to be an accelerated shelter development 
program. 

2. Tucson. Arizona 

Tucson has a residential population of 265,660, about the same as 
Orlando, Florida. The population density of 3000 people per square mile is 
below the 200 SMSA average. A look at the relative location of population and 
shelter indicates much the same posture as Orlando (Figure 9, Cumulative Distri¬ 
bution of Population and Shelter in Tucson, Arizona). 
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Percentage of Casualties Percentage of Population Sheltered 





Fallout Environment 



1 600 5.5 

2 2000 1.5 

3 3000 1 

4 5000 1 

5 10000 1 

6 30000 1 



Casualties Avoided by 
Assuming Improvable Spaces 
are Improved and Utilized. 


Excluding the SL Movement 
Restriction 


Fig. 8. Percentage of Casualties for Each Fallout Environment in Orlando. 
Florida (Total Population - 263,540) 
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300 



Radial Distance From Center of City (miles) 


Fig. 9. Cumulative Distribution of Population and Shelter 
is Tucson, Arizona 

According to NFSS Phase 1 data, 24 percent of the population can be shel¬ 
tered in spaces with protection factors ranging from 20 to 1000. Only a small 
number, 6 percent of the population, can be sheltered in spaces with 
PF )> 100 and the required area and volume. However, improving the existing 
spaces would provide shelter for an additional 16 percent with protection 
factor of at least 100. The cost of the improvements would be approximately 
$522,000, or $7.50 per space.. 

Figure 10, Percentage of Population Sheltered for Each PF Category in 
Tucson, Arizona, shows the shelter status of existing and improvable spaces. 
Although the fallout shelter status in Tucson is somewhat better than Orlando, 
the effect of movement alternatives is the same. (See Table V, Expected 
Casualties for Each Fallout Environment in Tucson. Arizona .) There just aren't 
enough shelter spaces available, even if improvement alternatives are incorpo¬ 
rated and spaces in PF Categories below 4 are accepted. 

Figure 11, Percentage of Casualties for Each Fallout Environment in 
Tucson, Arizona, again illustrates the reduction in casualties expected from 
assuming improvable spaces are improved and utilized. 
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Percentage of Population 
Sheltered in Improvable Spaces 



|~j Percentage of Population 

Sheltered in Existing Spaces 



PF Category 


Percentage of Population Sheltered for Each PF Category 
in Tucson, Arizona (Total Population = 265,660) 









Percentage of Casualties 


100 

90 

80 
70 

60 

50 
40 

30 

20 

10 

Fallout Environment 



Casualties Avoided by 
m Assuming Improvable Spaces 
are Improved and Utilized. 


NO. 


Fallout Environment 



1 600 5.5 

2 2000 1.5 

3 3000 1 

4 5000 1 

5 10000 1 

6 30000 1 


Excluding the SL Movement 
Restriction 


Fig. 11. lercentage of Casualties for Each Fallout Environment in 
Tucson, Arizona (Total Population - 265,660) 
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TABLE V 


I 

i 


Expected Casualties for Each Fallout 
Environment in Tucson, Arizona 


Fallout 1 

Environment 


Expected 

Casualties (%) 


Reference 

Arrival 

Movement 

Movement 

Constant 


Intensity 

Time 

Restricted 

Restricted 

Movement 

Movement 

(r/br) 

(hr) 

fo own SL 

to 2 Miles 

Speed 


Unrestricted 

600 

5.5 

89 

74 

74 


74 

2000 

1.5 

89 

74 

74 


74 

3000 

1 

94 

90 

90 

| 


90 

5000 

1 

97 

95 

95 


..95. 

10000 

1 

98 

97 

97 


97 

30000 

1 

99 

98 * 

98 


98 


3. Charleston, South Carolina 

This SMSA is somewhat atypical in that the population density—12,930 
people per square mile--is very large and the shelter population ratio of 0.33 
is rather low. The general trend is for high population densities to be 
associated with a relatively high shelter/population ratio (see individual 
points in Figure 3). Despite the combination of high population density and 
low shelter/population ratio, the relative location—with respect to the popu¬ 
lation center—of the population and shelter is essentially the same as in 
Orlando and Tucson. The shelter deficit exists throughoutjthe SMSA (Figure 12, 
Cumulative Distribution of Population and Shelter in Charleston, South Carolina) 
As might be expected, except for the SL restriction, there is no variation in 
the number of expected casualties regardless of the method used to allocate 
people to shelters. The expected casualties are shown in Table VI, Expected 
Casualties for Each Fallout Environment in Charleston, South Carolina and 
Figure 14, Percentage of Casualties for Each Fallout Environment in Charleston, 
South Carolina. 

The NFSS indicates that approximately 33 percent of the population can be 
sheltered in existing spaces (including Category 1). However, as in Orlando 
and Tucson, the great majority of spaces provide a protection factor of less 
than 100. Only 7 percent of the population in Charleston can be sheltered 
in spaces of 100 PF or better. Adding possible improvements would provide 
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Cumulative No. (Thousands) 


300 



Radial Distance From Center of City (miles) 


Fig. 12. Cumulative Distribution of Population and Shelter in 
Charleston, South Carolina ... 
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Fig. 13. Percentage of Population Sheltered for Each PF Category 

in Charleston, South Carolina (Total Population « 216,382) 
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Percentage of Casualties 





Casualties Avoided by 
Assuming Improvable Spaces 
are Improved and Utilized 



Fallout Environment 

No. 

i tr jrnm 

T A _£hii 

1 

600 

5.5 

2 

2000 

1.5 

3 

3000 

1 

4 

5000 

1 

5 

10000 

1 

6 

30000 

1 

* 




it 

Excluding the SL Movement 
Restriction 


Fig. 14. Percentage of Casualties for Each Fallout Environment 

in Charleston, South Carolina (Total Population * 216,382) 
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TABLE VI 


Expected Casualties for Each Fallout 
Environment in Charleston, South Carolina 


Fallout Environment 

Expected Casualties (7.) 

Reference 
Intensity 
. (r/hr) 

Arrival 

Time 

(hr) 

Movement 
Restricted 
to own SL 

Movement 
Restricted 
to 2 Miles 

Constant 

Movement 

Speed 

Movement 

Unrestricted 

600 

5.5 

83 

70 

70 

70 

2000 

1.5 

83 

70 

70 

70 

3000 

1 

92 

89 

89 

89 

5000 

1 

97 

95 

95 

95 

10000 

1 

99 

98 

98 

98 

30000 

1 

99 

99 

99 

99 

shelter 

for another 

6 percent of 

the'population. 

The total cost 

of improve- 


ments would be $130,000 or $8.44 per space. 

4. Mobile. Alabama 

Mobile has a residential population of 273,942, which is the second largest 
of the eight SMSA's being considered. The population density is the lowest at 
1320 people per square mile. Although the SMSA has an overall shelter deficit 
of 54 percent, it is different from the previous three SMSA's in that a shelter 
surplus does exist within a two-mile radius of the population center (Figure 15, 
Cumulative Distribution of Population and Shelter in Mobile, Alabama). In 
fact, the area outside the two-mile radius has shelter for only 8 percent of 
the population, and 62 percent of the total population is located in the same 
region. Therefore, about 6 percent of the population must ' -s moved "downtown" 
to the available shelter spaces. The expected numbers of casualties for the 
various movement restrictions are shown in Table VII, Expected Casualties for 
Each Fallout Environment in Mobile. Alabama . As can be seen, there is little 
variation in casualties according to the method of allocating people to shelter, 
with the exception of the extreme case where the population is restricted to its 
own standard location. 

Using spaces located by the NFSS, 48 percent of the population can be 
sheltered in spaces with PF ranging from 20 to 1000. (Figure 16, Percentage 
of Population Sheltered for Each PF Category in Mobile, Alabama). Approximately 














0 


Total Population ■ 273,942 











Percentage of Population Sheltered 



Fig. 16. Percentage of Population for Each PF Category in Mobile 
Alabama (Total Population - 273,942) * 
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TABLE VII 


Expected Casualties for Each Fallout 
Environment in Mobile, Alabama 


Fallout 

environment 

Expected Casualties (7.) 

Reference 

Arrival 

Movement 

Movement 

Constant 


Intensity 

Time 

Restricted 

Restricted 

Movement 

Mbvemen t 

(r/hr) 

(hr) 

to own SL 

to 2 Miles 

Speed 

Unrestricted 

600 

5.5 

78 

.._ . 54 _ 

—. 54 

54 

2000 

1.5 

78 

54 

54 

54 

3000 

1 

91 

83 

83 

83 

5000 

1 

96 

. 93 

93 

93 

10000 

1 

97 

96 

96 

__ 96 _ 

30000 

i 

99 

99 

99 

99 


15 percent of the population can be sheltered in spaces having FF of at least 
100. Improving existing spaces would provide shelter for an additional 9 per¬ 
cent of the population at a cost of $170,000 or $6.90 per space. 

5. Macon. Georgia 

Macon has very nearly the same shelter/population ratio, 0.50, as Mobile; 
however, the population density is much higher. The total population is 141,249, 
and the population density is 4650 people per square mile. Macon, like Mobile, 
has a shelter surplus near the population center of the SMSA and very few 
spaces outside a two-mile radius of the center (Figure 18, Cumulative Distribu¬ 
tion of Population and Shelter in Macon, Georgia). 

According to NFSS data, 50 percent of the population can be sheltered, and 
20 percent of the population would have at least a PF of 100. An additional 
14 percent of the population could be sheltered in improvable spaces. Ihe cost 
of the improvements would be $94,000 or $4.75 per space. 

The distribution of shelter spaces by PF category shows an increase in the 
fraction of spaces having a PF of 100 or better over SMSA's previously described 
(Figure 19, Percentage of Population Sheltered for Each PF Category in Macon, 
Georgia). As was the case for all previously examined SMSA's, the largest 
increase in casualties, for all movement restrictions—as shown in Table VIII, 
Expected Casualties for Each Fallout Environment in Macon. Georgla --occurs when 
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Radial Distance From Center of City (miles) 

Fig. 18. Cumulative Distribution of Population and Shelter in 
Macon, Georgia 









Percentage of Population 
Sheltered in Improvable Spaces 



PF Category 

Fig. 19. Percentage of Population Sheltered for Each PF Category 
in Macon, Georgia (Total Population - 141,249) 



Percentage of Casualties 
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1 2 3 4 5 6 


Fallout Environment 


Casualties Avoided by 


Assuming Improvable Spaces 
are Improved and Utilized. 



Fallout Environment 

No. 


4-ftr). 

1 

600 

5.5 

2 

2000 

1.5 

3 

3000 

1 

4 

5000 

1 

5 

10000 

1 

6 

30000 

1 

Excluding the SL 
Restriction 

Movement 


"k 

Fig. 20. Percentage of Casualties for Each Fallout Environment 
in Macon, Georgia (Total Population ■ 141,249) 


TABLE VIII 


EXP 


ected Casualties for Each Fallout 
Environment In Macon, Georgia 


Fallout Environment 

Expected Casualties (7°) 

Reference 

Arrival 

Movement 


Movement 

Constant 


Intensity 

Time 

Restricted 


Restricted 

Movement 

Movement 

(r/hr) 

(hr) 

to own SL 


to 2 Miles 

Speed 

Unrestricted 

600 

5.5 

77 


46 

46 

46 

2000 

1.5 

77 


46 

46 

46 

3000 

1 

85 


73 

73 

73 

5000 

1 

87 


84 

84 

84 

10000 

1 

90 


86 

86 

86 

30000 

1 

97 


97 

97 

97 

the fallout reference 

intensity changes from 2000 

r/hi to 3000 r 

/hr. The reason 

for the sharp increas 

e is due to 

i 

the 

ineffectiveness of Category 

1 spaces at 

3000 r/hr. 

Of the to 

cal number of spaces available, more than half of them are 

Category 1 

spaces. A 

characteristic 

i 

of Macon not 

encountered in 

the previously 

considered 

SMSA's Is 

the increase 

in 

casualties at 

the 30,000 r/hr fallout 

reference 

intensity. 

The SMSA's, 

previously considered, experienced rather 

constant marginal increases over 

the 

range from 3000 r/hr to 30, 

000 r/hr. 

Macon, however, has a 

relatively large number of Category 5 and 

6 spaces (enough 


for 12 percent of the'population) which are ineffective at reference intensity 
of 30,000 r/hr, thereby causing the increase in casualties when I q » 30,000 r/hr 

6. Pittsfield. Massachusetts 

Pittsfield has the smallest population in the sample (73,839), and the 
second smallest population density (1420 people/square mile); however, its 
shelter/population ratio is the third largest (0.66). A shelter surplus exists 
within a one-mile radius of the population center (Figure 21, Cumulative Distri¬ 
bution of Population And Shelter in Pittsfield, Massachusetts). It is also 
the only SMSA with enough improvable spaces to erase the shelter deficit. By 
using the improvable spaces, 47 percent of the population could be sheltered 
in spaces having a PF of at least 100. (Figure 22, Percentage of Population 
Sheltered for Each PF Category in Pittsfield, Massachusetts), 




















lative No. (Thousands) 



Radial Distance From Center of City (miles) 

Fig. 21. Cumulative Distribution of Population and 
Shelter in Pittsfield, Massachusetts 
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Percentage of Casualties 



Casualties Avoided by 
Assuming Improvable Space 
are Improved and Utilized 


No. 

Fallout Environment 

1 

600 

5.5 

2 

2000 , 

1.5 

3 

3000 

1 

4 

5000 

1 

5 

10000 

1 

6 

30000 

1 


Excluding the SL Movement 
Restriction 


Fallout Environment 

Fig. 23. Percentage of Casualties* for Each Fallout Environment 

in Pittsfield, Massachusetts (Total Population ■ 73,839) 
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As in all the other SMSA's, intense fallout environments produce a high 
casualty rate as shown in Table IX. 


TABLE IX 

Expected Casualties for Each Fallout 
Environment in Pittsfield. Mass. 


Fallout 

Environment 

Expected Casualties (%) 

Reference 

Arrival 

Movement 

Movement 

Constant 


Intensity 

Time 

Restricted 

Restricted 

Movement 

Movement 

(r/hr) 

(hr) 

to own SL 

Co 2 Miles 

Speed 

Unrestricted 

600 

5.5 

64 

43 

43 

43 

2000 

1.5 

64 

43 

43 

43 

3000 

1 

72 

66 

66 

66 

5000 

1 

81 

81 

81 

81 

10000 

1 

87 

87 

87 

87 

30000 

1 

. 90 

90 

90 

90 


The variation in expected casualties is quite large when comparing condi¬ 
tions where some movement is allowed and those conditions where movement is 
restricted to the SL. For example, in the fallout environment which would be 
exceeded in only about 15-30 percent of the attacks (i.e., I q ■ 2000 r/hr, 

T^ « 1.5 hrs), approximately 64 percent of the population become casualties 
if movement is restricted to within the SL while only 43 percent are casualties 
with any other mo\ement restriction. The point being made is that casualty 
computations are sensitive to movement restrictions in the majority of expected 
"fallout environments. 


Thus far the sensitivity has been evident only when movement is restricted 
to the SL compared to movement outside the jSL. In the remaining SMSA's, 
casualties becomt sensitive to all four movement restrictions. 


7. Waterloo. Iowa 

Waterloo is characterized by a shelter surplus over the entire SMSA 
(Figure 24, Cumulative Distribution of Population t.nd Shelter in Waterloo, Iowa). 
The characteristics of Waterloo are: population--122,482; population density-- 
2120; shelter/population ratio--1.10. Being fortunate enough to have enough 
spaces, hovever, also creates new problems. In all the previously discussed 
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100 



1 2 3 4 5 6 


Fallout Environment 

D No Movement Restriction, No Movement in Fallout. 

® Casualties Added—Movement Restricted to 2 MPH, Movement in Fallout 
is Allowed. 

■ Casualties Added—Movement Restricted to 2 Miles, No Movement in 
Fallout. 


S Casualties Added—Movement Restricted to Own SL, No Movement in Fallout. 


Fig. 25. Sensitivity of Casualty Calculations to Movement Restrictions— 
Waterloo, Iowa (Total Population - 122,482) 
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30000 1 
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Algorithm Solution 


Fig. 27. 


Percentage of 
in Waterloo, 


Casualties for Each Fallout Environment 
Iowa (Total Population - 122,482) 





SMSA's, Che shelter deficit was such that there was no serious problem in 
utilising all the spaces. In fact, the available shelter could be treated as 
if the population has unrestricted movement within the SMSA. In Waterloo, with 
its shelter surplus, the expected number of casualties is quite sensitive to 
any change in movement restrictions. For example, consider the results of the 
allocation model as shown in Table X. Only in the most severe fallout environ¬ 
ments are the casualty calculations insensitive to the movement restrictions. 


TABLE X 


Expected Casualties for Each Fallout 
Environment in Waterloo. Iowa 


Fallout 

environment 

Expected Casualties (7.) 

Reference 

Arrival 

Movement 

Movement 

Constant 


Intensity 

Time 

Restricted 

Restricted 

Movement 

Movement 

(r/hr) 

(hr) 

to own SL 

to 2 Miles 

Speed 

Unrestricted 

600 

5.5 

83 

49 

0 

0 

2000 

1.5 

83 

49 

25 

0 

3000 

1 

87 

54 

42 

22 

5000 

1 

89 

64 

58 

57 

10000 

1 

91 

68 

68 

68 

30000 

1 

96 

83 

83 

83 


These same results are shown graphically in Figure 25, Sensitivity of 
Casualty Calculations to Movement Restrictions--V'aterloo, Iowa. Quite obviously, 
some movement is necessary in order to utilize shelters and avoid large numbers 
of casualties. As will be shown later (in Section V and Figure 31), because 
of the gradual buildup of fallout particles, some movement in a fallout field 

is permissible for reference intensities as high as 5000 r/hr,—Restricting- 

movement seriously degrades the effectiveness of existing shelter spaces. 

8. Richmond. Virginia 

This SMSA has the largest population (408,494) of all of the SMSA’s in the 
Sample. It also has a 60 percent shelter surplus, (Figure 28, Cumulative Distri¬ 
bution of Population and Shelter in Richmond, Virginia); and by incorporating 
shelter improvements, approximately 60 percent of the population could be shel¬ 
tered in spaces having PF of at least 100. Here, as in Waterloo, casualty 
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Cumulative No. (Thousands) 



Radial Distance From Center of City (miles) 


Fig. 28. Cumulative Distribution of Population and Shelter in 
Richmond, Virginia 
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it 

Fig. 30. Percentage of Casualties for Each Fallout Environment 
in Richmond, Virginia (Total Population * 408,494) 
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calculations are sensitive to proper shelter utilisation. For example, consider 

a fallout environment, I ■ 3000 r/hr, T. * 1 hour. If movement is. restricted 

o A 

to within standard locations, the expected casualties will be 81 percent. How¬ 
ever, if movement is unrestricted (i.e., assume strategic warning) there will 
be one percent casualties. Obviously, these are extreme examples. Using an 
optimum allocation plan, such that casualties are minimized subject to a move¬ 
ment speed restriction, would result in approximately 38 percent casualties. 
Using the two-mile restriction, which is similar to the optimal plan except 
that it allows no movement in fallout, would result in 68 percent casualties-- 
a 30 percent increase. The significance of limited movement in fallout is quite 
evident. The expected casualties at various other fallout environments are 
shown in Table XI. 



TABLE XI 


Expected Casualties for Each Fallout 
Environment in Richmond, Virginia 


Fallout Environment 


Expected 

Casualties (%) 


Reference 

Arrival 

Movement 

Movement 

Constant 


Intensity 

Time 

Restricted 

Restricted 

Movement 

Movement 

(r/hr) 

(hr) 

to own SL 

to 2 Miles 

Speed 

Unrestricted 

600 

5.5 

66 

47 

0 

0 

2000 

1.5 

66 

47 

36 

0 

3000 

1 

81 

68 

38 

1 

5000 

1 

88 

81 

50 

45 

10000 

1 

90 

88 

87 

70 

30000 • 

1 

93 

91 

90 

90 



V. DATA ANALYSIS 


V. The Relationship Between Time in Fallout and Protection Factors 

An analysis of ERD shows interesting trade-offs between time in fallout and 
:he protection factor of the shelter that is entered. Even when the reference in- 
;ensity is 5000 r/hr and the time of fallout arrival is only one hour, considerable 
:ime in fallout can be compensated for by a relatively small increase in the protec- 
:ion factor of the shelter eventually entered. For example, a person traveling 15 
ulnutes in such a fallout field (at 2 mph) and reaching a fallout shelter with a PF 





4 


of 91 would sustain no more dose than a person who was in a shelter with a PF of 
70 at the time fallout commenced. Of course, this type of comparison considers, only 
gamma radiation. Other types of radiation might also have to be taken into considera¬ 
tion, but their effect should be negligible if inhalation of particles is avoided. 

Figure 31, Trade-Off Between PF and Time in Fallout, shows the trade-offs of 
PF and time in fallout for the cases where the reference intensities are 5000 r/hr, 
_3000 r/hr ^ and 2000 r/hr. The latter has an associated time of arrival of 1.5 hours; 
the others have a time of arrival of one hour. The curves are based on a buildup 
time of twice the time of arrival and a linear rate of intensity buildup. 

In the case of reference intensity of 10,000 r/hr, there is little opportunity 
for trade-off between PF and travel in fallout. A person 

environment would receive a dose greater than 200r ERD if 

% 

as ten minutes, regardless of the type of shelter he even 


unsheltered in such an 
h^e were outside as much 
ually entered. 


Relationship Between Reference Intensity, Protection Factor, and Maximum ERD 


The ERD calculations also provide a means for comparing reference intensity, 
protection factor, and Maximum ERD. Figure 32, Trade-Offs Between PF and Fallout 
Environment, shows these relationships for the reference intensities considered 
in this volume. The dashed lines show the effect of times of arrival of more than 
one hour. 

' » 


VI. CONCLUSIONS AND RECOMMENDATIONS 


A. Discussion j 

Shelter utilization analyses necessarily include assumptions concerning 
variables that are difficult to measure. These variables affect utilization. Among 
the more significant variables are: ... 

1. Knowledge on the part of people about the shelter which they should use 
(shelter assignment plan, discipline, crowd control, etc.). 

2. Time required to reach shelter (warning time, start-up time, travel time, 
searching time, interference of traffic flow, etc.). 

3. Distance to shelter (influences both time and knowledge, for the knowledge 
may be less if distances are long). 

In a complete analysis, relations among these three and other variables must 
be derived (for they are obviously not independent), and these in turn related to 
the shelter utilization pattern. 
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In the analyses reported here, knowledge, time, and distance to the shelter 
were introduced in a semi-quanti^ative way by varying "movement restrictions" as 
follows: 

1. Movement Restricted to Own SL 

Knowledge very great, and time and distance allowance great enough to 
allow best shelter utilization within one's own SL. 

2. Movement Restricted to 2 Miles From Own SL 

Knowledge very great, and time and distance allowance great enough to 
allow best utilization of any shelter within one's Standard Location (SL), or 
a SL not more than two miles away (rectangular distance). 

3. Movement Unrestricted Within SMSA 

Knowledge very great, and time allowance great enough Lu allow best shel¬ 
ter utilization, at long distances anywhere in the SMSA. 

4. Optimum Shelter Utilization (Transportation Algorithm) 

Knowledge exact, to allow optimal shelter utilization (minimum casualties) 
anywhere within the SMSA, including a prior calculation of dose received during 
transit in fallout and in the primary shelter, all before cr during fallout 
arrival. Speed was restricted to 2 tnph. 

We have a mathematically optimum solution in Movement Restriction 4. A similar 
ihelter utilization pattern could be generated by other assumptions, operationally 
tore realistic--as seen In Table XII, Summary of Casualty Estimated (% of Population) 
'or Reference Intensity of 2000 r/hr and Time of Arrival of 1,5 Hours --for unrestricted 
lovement, for example. For more severe environments, the two casualty figures (hence 
helter utilization) become identical. This case and unrestricted movement offer a 
■enchmark of ideality against which to measure other shelter utilization patterns. 

. Major Conclusions 

Efficient shelter utilization appears to be a critical problem when the shelter 
paces exceed population; inefficient utilization can lead to a large number of 
voidable casualties. This was concluded from the observation of the rapid increase 
n casualties for a specified attack environment as one restricts movement progres- 
ively from the SMSA down to the SL in which the shelteree finds himself at the time 
f warning. For the most shelter-abundant cities (Waterloo and Richmond in Table 
II), the percentage of casualties under a 2000 r/hr fallouc environment is 
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progressively reduced from a tvo-city average ranging between 757. and 307. for limited 
movement (SL and Constant Movement Restriction in Table XII) to zero for movement 
unrestricted within the SMSA. 


TABLE XII 

Summary of Casualty Estimates (% of Population) for 
Reference Intensity of 2000 r/hr and Time of Arrival of 1.5 Hours 


SMSA 

Shelter/ 
Population 
Ratio (%); 
Categories 

1 and 
Better 

Movement 

Restricted 

To Own SL; 

Movement 

Completed 

Before 

Fallout 

Constant Movement 
Speed of 2 mph 
Restriction 
Movement Permitted 
in Fallout (Trans¬ 
portation 
Algorithm) 

Movement 
Unrestricted 
Within SMSA; 
Movement 
Completed 
Before 
Fallout 

Orlando 

20 

94 

78 

78 

Tucson 

24 

89 

74 

74 

Charleston 

33 

83 

70 

70 

Mobile 

48 

78 

54 

54 

Macon 

50 

77 

46 

46 

Pittsfield 

66 

64 

43 

43 

Waterloo 

110 

83 

25 

0 

Richmond 

160 

66 

36 

0 


For the most shelter-poor cities (Orlando and Tucson), the corresponding reduc¬ 
tion is from a two-city average of 917. casualties to a two-city average of 767. casu¬ 
alties, a feasible savings of 15% of the population who would otherwise be casualties 
This is to be compared to the 757. savings to the more shelter-abundant cities. The 
lack of sensitivity to movement restrictions in shelter-poor cities generally is due 
Co the fact that there are enough people in any standard location to completely 
utilize available shelter in that SL. Thus, even the most restrictive shelter utili¬ 
zation constraint leads to good utilization. 

Within broad limits, tne higher the ratio of shelter is to population, the 
greater is the payoff in casualties saved by shelter utilization planning. 

Insofar as movement distance restrictions are derived from time available for 
sheltering, one can deduce the substantial payoff potential in effective warning 
means, shelter planning, and indoctrination of the population--particularly in 
shelter-abundant areas. 

C. Secondary Conclusions 

1. Movement in fallout to a better shelter theoretically is warranted in 

many circumstances. The quantitative trade-offs between movement time in 
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fallout and protection factor of the shelter are shown in Figure 31. Un¬ 
feasibility of applying this principle in the transattack environment is 
questionable, however. It is a special case of remedial movement, but with a 
very short decision time and a changing environment. 

2. Neither city size nor population density (one of the city selection' criteria) 
which might be supposed to have an effect on juxtaposition of people arid shelter, 
could be shown to be significant in this analysis. This is concluded from 
examining the difference in casualties in each attack environment, between 
Movement Restriction 4, (movement in fallout, but limited movement speed), and 
Restriction 3 (complete freedom) as a function of city size and population 
density. Small effects which likely exist either are swamped by the more 
critical shelter/population ratio, or are undetectable from the "noise" of 
variable PF distributions in cities otherwise similar (size and population 
density). The casualties differ between Restrictions 3 and 4 only for Waterloo 
(low density, small size) and Richmond (high density, large size) at the 
intermediate attack environments. 

I). Recommcndat ions 

This analysis has shown the critical sensitivity of casualty estimates to 
shelter utilization, and therefore, to the corresponding movementrto-shelter 
iissumptions, particularly for shelter-abundant areas. Accordingly, further research 
is recommended to give a more quantitative insight into the problem of movement-to- 
shelter (time, training and discipline, and distance to move). Better understanding 
can lead to' more effective operational plans and programs in terms of preattack 
shelter assignment planning, warning requirements, training, etc. 
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Appendix A 

Development of Probabilistic Fallout Environments 

1. INTRODUCTION 

A. Cumulative Distributions Mathematically Defined 

Given a simulated nuclear attack, the pattern of the attack is used to compute, 
on a probability basis, the severity of the fallout environment at each of 200 SMSA's 
throughout the U. S. One-hundred attack patterns comprise one RISK attack. Thus for 
one RISK attack, each SMSA will be subjected to 100 possible fallout environments. 

Let I and T^ dnote the Reference Intensity and Time of Arrival, respectively. 
Each component is divided into subgroups: 

A similar distribution is derived for time of arrival T.: 


i 

I 

0 

t a 


(r/hr) 

(hr) 

1 

0-100 

0-1 

2 

101-300 

1-2 

3 

301-1000 

2-3 

A 
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5 

3001-10000 
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6 

10001-30000 
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If X. . denotes the number pf times I occurs in the interval for the i 1 * 1 

ij 0 th 

SMSA then the probability of occurrence for the 1 interval is estimated by 
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where f represents the number of attack patterns. Similarly, if denotes the 
number of occurrences of in the i interval for the SMSA, we have 
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The discrete cumulative distribution function for each variable is developed by 
summing the individual probabilities. 
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Tin.' Populat ion to Shelter Algorithm 


I 

I. INTRODUCTION 

■I 

A. Expression of Content . /'--- •.... 

This appendix describes in detail the model used to estimate casualties within 
an SMSA. Also included is the computer program used in calculating the Maximum ERD 
associated with movement between standard locations. 

1. Development of a Mathematical Model to Calculate Maximum Equivalent 
Residual Dose (Stage 1) 

The Model and Program described in this section is a modified version of 
the Mainline Model and Program discussed in Appendix of Volume III, "A 

Generalized Sensitivity Analysis of CD Systems." The greatest difference between 
this algorithm and the Mainline Model of Appendix B of Volume III is in the input/ 
output part of the program. By necessity each being oriented to a different re¬ 
search objective. 

The equation describing ERD at t hours is the sum of the reparable and 
irreparable dose over finite time increments. 



Fig. B-l. Illustration of Finite Different Formulation 
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T. *nd T denote the time of arrival and time of cessation of fallout 

Ac 

respectively. 

* Hr 

Let Dp represent the reparable dose, and the irreparable dose for the 
p** 1 interval. Then . 


D* - [1 - (a) (b) ] D* l + (1-f) ^11^- 


(B-l) 


** ** 

D„ - D n + (f)(r)(a) 
P P"1 


(B-2) 


where 


(1) a is the constant time increment in hours. 

(2) b is the recovery rate (0.1%/hr.). 

(3) f is the nonreparable part of the dose. 

(4) PF is the protection factor for the shelter used. 


r (t - a/2 - T a ) 


T. < t < T 
A — — c 


(B-3) 


T = 2T. 
c A 


r ■ a t 


, t > T 

' f' 


(B-4) 


The ERD for the p time increment is 


* ** 

D - (D + D ) I , 
P P P o 


(B-5) 


and the Maximum ERD occurs when D > D ., 

P P+1 


Ihe previous derivation is based on the assumption that a person enters a 
shelter with protection PF before fallout arrives; i.e., if the time required 
to move from standard location i to j is T^^, then the assumption is that 
T^j < T a . If T fi > T a , then in Equation B-l, PF must equal 1 for T { , > t . 


l ij A 


When T^ becomes less than or equal to t , (T^ < t) the appropriate PF is 
empLoyed and the original derivation is valid. 


The computer program for Stage 1 is described in Sections 3, 4, 5 of this 
Appendix. 
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Th«A« s i z —se-u x, .x!« \ (St aze 2) 


The assignment n,>»e! i« t-a«e 1 on a linetr ft •,.raining technique called 
*.' , h i • » M* * - * . " Itiis technique is most commonly employed in 

'.it ' t* * l:.e cost of shipping a homogeneous product manufactured in "n" mill 
4 *\ch in a different geographical location, and shipped to consumers at "m" 
different destinations. The analogy is made that people represent the product 
located in "n'Lstandard locations and fallout shelter spaces represent the 
warehouse capacities at "m" standard locations. The cost variable is the 
equivalent residual dose received if one moves from standard location i to 
standard location j where j has a protection factor of k . ^ 

- Symbolically, the model formulation is: ." - *.I------ r 


Let represent the ERD defined previously: |'|J' 

s ijk ■ f<R - T *- v - s ij- pr k> < b -m 

R "» fallout reference intensity (r/hr) 

T. — arrival time of fallout (hr) 

A 

V ~ movement speed (mi/hr) 

distance from SL i to SL j (mi) 

PF k ~ protection factor. 

■ *' 

The model which calculates s,„ is defined in Section A.l. Then assume the 

ijk 

ERD obtained in moving from any SL i (i ■ 1, ..., n) to any SL j (j*>l, .... 

✓ 

with a Drotection factor category k (k ■ 1, 8) Is known. Thus, we have 

a 3-dimensional matrix (n x m x 8 ) with each element representing a."cost" 

t Hi >' . 

[or Maximum ERD. To simplify the matrix, let jk represent the j destination 
th' 

and k protection category. What this means is that each protection category 
is treated as a separate and distinct destination, where 


jk - 8(J-1) + k . 


(B-7) 


For example, let the cost matrix be represented by D and let s represent 

the Maximum ERD for a person leaving standard location 3, moving to standard 
location 5, and entering a protection factor category 8 . The element in D 
corresponding to s 3 5 g is where 


i - 3 

jk - 8(5-1) + 8 - 40 . 


(B- 8 ) 











There fere 


*1 t c — » 1 A * nd S , - D . 

3,f ,8 sr 3 . -*0 3,-0 

To summarize briefly the derivation of the cost matrix: 

(a) The maximum ERD is calculated in Stage 1 of the model for each 
element of a 3-dimensional matrix where s ^ represents each element 

—in the matrix. ;; ~ •' .“""",7; . .. 

(b) In order to fit the 2-dimensional cost matrix scheme of the transporta¬ 
tion algorithm, a linear transformation was performed on each element 

whereby s. became s,,,. v . 

ijk Kjk) 


With the cost matrix in the desired form, the constraint equations and 
objective function may be defined: 

(a) The population allocated from origin i to destination jk(X . .) 

* \ J K/ 

must not exceed the initial population at the origin, A^. 


Z x (jk) < A , (i « 1.n) 

jk-1 1 1 


(B-9) 


X i(jk) - 0 * 

(b) The population sheltered at destination jk must not exceed the shelter 
capacity at the destination, B. ... 


X i( jk) - B ( jk) * Ok “ 1. .... 8m) 


(B-10) 


X i(jk) - 0 * 


(c) The objective function to be minimized is the number of casualties, 
subject to the previously defined population and shelter space con¬ 


straints. 


n 8m 

min z « Z Z c..,.. x. .. 

i-1 jk-1 l ^ k) ^ k) 


i ( jk.) 


°’ d i(jk) - 200r 
l * d i(jk) > 200r 


(B-ll) 
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ISL Number of standard ! >c a t i s in the 'jMSA 

LMAX Maximum number of PF categories 

Z Radiation decay constant 

R1 Reference intensity 

TA Arrival tine of fallout 

VEL Movement speed 

B Recovery rate 

F Recoverable fraction 

DMAX Maximum acceptable ERD 

PF1 Primary PF 

PF2 Secondary PF 

PF3 Equivalent PF in third shelter period 

PF4 Equivalent PF in fourth shelter period 

T1 Time in primary shelter 

T2 Time in secondary shelter 

T3 Time in third shelter 

III SMSA code 

IP Number of populated standard locations 

IS Number of standard locations with shelter spaces 

IC IS x IP 

RMAX Peak reference intensity of buildup function 

D1 Reparable portion of total dose 

D2 Irreparable portion of total dose 

DD Equivalent Residual dose (ERD) 

DP Maximum ERD 

TMM Cumulative time since time of arrival 

L. . .. PF category. 1 ~ ~ “ 
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I BEGIN LOOP 
TO CALCULATE 
MAX. ERD FOR 
MOVEMENT FROM 
3l(i) TO SL(j) 


CALC. 

<Jl) DISTANCE 
TRAVELED 
(2) TIME IN 


Exit 
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Fortran Listing 


76655 ISS9XI 1C56 

• . . ..SC3NCV SiM . .RH-QCD 

• xe r; 

SECTION I - CALCULATION OF MO VEMFn T TIMES 

OIMENSICN XI100» »VI ICC) v n(l01*01(10 I*02(101,CPIlC>,TMtlCI 
DIMENSION PFiUO>,PF2(iOl ,S110CC1,JRSAC 1100) 

DIMENSION, DDPIAG.AC.B),TFRm(AG,AC) ,PflOO) 

rewinc g 

IGNCRF =C 

1 CONTINUE 
11=1 

RFAC INPUT TAPE 5,3C1,I$L,LMAX,? 

IFUSl 171A.71A.2 

2 READ INPUT TAPE 5,5C1,R1,TA.VEL 
READ INPUT TAPE 5,5C2,B,F,DMAX 

1C CONTINUE 

READ INPUT TAPE 5 , A , PF1(II,PF11 2 I,PFIt 3 I,PF1(A ) ,PF1(5» 

READ INPUT TAPE 5,A,PF1{6) ,PF1C7),PF1I8) ,PF1(9)«PF1UC» 

RE AC INPUT TAPE 5»A»PF2I1),PF2C2)»PF2I3)»PF2(A),PF2(5I 
RF AC INPUT TAPE 5 , A , PF21 6 I,PF2C7),PF2 I A),PF219),PF211C» 
READ INPUT TAPE 5 , A, PF 3, PF A , T1, T2, T 3 
3C1 FORMATI2I5*F1C.2) 

5 C 1 FORMAT I3F1C.2) 

5C2 FORMAT I 3F1C.5) 

A FORMAT (5F1C.2) 

WHITE CUTPUT TAPE 6,8C 1 ,R1.TA,VFL,PF1f1) 

8C1 FORMAT (AF12.2) 

RE AC INPUT TAPE 5,7C7 

707 FORMAT(ACHI ) 

WRITE OUTPUT TAPE 9,7C7 
WRITE CUTOUT TAPE 6,707 
READ INPUT TAPE 5,7C5.III,IP,IS,10 
WRITE CUTPUT T AP r 6,7C6,1P , IS,IC 
WRITE CUTPUT TAPE 9,7C6,1P , IS,IC 
7C 5 FORMAT! 17,1 13,11C, 11G) 

?C6 FC R_M A_T ( 7HC 2 1111 3 I 1 3 ,110, I 1C 1 
IF( IGNORE)7061,7061.7C62 
706 1 DC 703 M* 1, ISL 

REAC INPUT TAPE 5,701,JRSAC(MJ,X{M),YIM} 

701 FORMAT!18,F12.C.F10.0) ’ 

7C3 CONTINUE 
706 2 CO NTINUE 

RMAX * (2.0*T A !•*{-2) 

7C8_C0NTINUE . 

WRITE CUTPUT TAPE 6,8CC 

flOC_fCRMA_T (5 5H1 I J PF JJ PRCB MAX ERD 

DO 305 1*1,ISL 

PC .305 J* I., ISO _ .. . . _ 

5C9 DO * C 
Dll * C 
D22 * C 
TMW * C „ 

DPP * C 

_t * ($cRTFnxm -x uj i”2) ^Rif mm-vij i Lt-2ii/m_. ... 

T * T»0,062 
§90 CONTINUE 
5CA CONTINUE 
5t£ .CONTINUE ... 


TIME) 



1 


T « T-TA 

SECT ICN I - CALCULATION OP MOVEMENT TIMES 

JF_CT1 5C_7» 5C81_508 

5C7 T * 0. 

_5C8 CONTINUE_ 

IF(T) 599,599,510 

510 TMM » TMM»Q, 1 _ 

RM *(RMAXMTMM-C.C5)/TA')*C.1 
IF{{RM» 10. l-RMAX>515*515,514 
514 l.-Z) 

__ _A2 a (T M M- . 1 }••{ l.-Z) 

RM«(A1-A2)/Il.-Z) 

-515 011_» y_l- .0-0.1»B 1_*C11 )♦ I (1.0-F >HI 

022 « 022*IF *RM) 

_CC_ = (011*022)»Rl 

IFIDD-CPP) 513,511,511 . 

_5ii opp = cc 

IF!DPP-DMAX)513,599,599 
513 IF I TMM -(T-0.05)1 51C.599.599 
599 CONTINUE 

_C _ SECT IC N_H I 

5 CONTINUE 
_7 _ OC 101 L = l.LMAX 
OIL I * CD 

_ DHL) » Dll _ 

02 IL) '= 02 2 
_DP(L) * OPP 
TM(L1 * TMM 

_ IFICP(Ll-OMAX) 19,999,999 

15 IF(TMM-T A) 20,30,30 

_ TM{l)*C.5_ 

RML =(RMAX*'( ITMll )-C.*25j/TA) )«0.5‘ 

Ol(L) * (l.C-B/2.C)»01(U + (l.C-F)*RML/PFl(U 
02 IL j = D2(L)*F»RML/PF1(L) 

OILI = R1*IC1(L)*02 t L ) 1 
IF (OIL )-DP (D) 22,22,21 

21 TERM(I ,J) = TMIL) _ _ 

OP (L) * DILI 

IF I CP I U-700. 122,999,999 
22 TAL * TA-0.25 I 

IF|TM(U-TAU2C,3C,30 
30 TMILI = TMIL) +6.0 

_ M *_TMIL ) •• I 1 .0—Z1 _ _ _ < 

A2 s (fMIL)-6.C)*•Il.C-Z) 

RML * (A1-A21/Il.C-Z) 

OKU = (l.C-6.*B)*01 C L I ♦ C l.O-F) •RML/PF1IL ) 

. 02 IL ) = 02(U*F*RML/PFl(U 

OIL) * R1MD1 (LKC2IL ) ) 

_ LLl-SlU-DfJU I,32i 32.31. . 

31 TERM I I,J) = TM(L ) 

_ OP(L)_ = 0(L I _ 

IF (CPU )-700.)32,995,999 

_ U IF(TM(L) -T1)3 C ,40,40 

40 CONTINUE 

4C1 TM(L) » TM (L ) ■♦■24. C _ 

A1 * TMIL)*»<l.C-Z) 

_1 -6.Q)** LUfl-i ) . . 

RML - ( A1-A21/(1.C-Z) 

I___ S££UflH_I ~ CALCULATION DF MOVEMENT TIN1S 





. . . 01 CL) * (l.£-2*.0»6)*01(M 

01(L) = 01 (L I ♦ U l.C-F )*RML >/PF?(U 
02(U * D2(L)*F*RML/PF2d) 

0 ( 1 ) = Rl*( 0 l(L)*C 2 d) ) 

IP (OIL )-OP(l) )A2,A2,Al 
Al TERM(I,J)=TM(L) 

OP(L) = 0(L ) 

IF(CP(L)-700. >A2,99<9,991 
A2 TAL = TUT2 

IF(TM(LJ-TAl) ACi.A9.A9 
A9 IF(T3)999,999,50 
5C T*M L ) = TM(L )*168.0 
A1 = TM(t)»*( l.C-Z) 

A2 = (TMd)-6.0)*»(l.C-Z) 

RVL = (A1-A2)/(X.C-Z) 

Old) = ( 1 .C- 16R.C-B )*01 { L )♦ ( l.C-F ) *RMl/PF3 
02(1) = 02(L)♦F*RML/PF3 
0 ( L ) = R1»(01(L)*C2(L) ) 

IF(C(L )-OP(U ) 52,52,51 
51 TERM(1,J j = TP{L) 

OP(l) = D(L) 

IF(CP(L)-700.)52,999,999 
52 TAL = T1+T2+T3 

IF(TM(t)-TAL)5C»39,59 
5 9 CCnTINUF 

6Q TM(l) = TM(U*720. 

A1 * TM(L)*«(l.C-Z) 

A2 * (TM(L)-6.0)**(1.C-Z) 

Rml = (A1-A2)/(l•C-Z) 

Ol(L) = ( 1.0-720. C»B)*C1(UM1.C-F)*R'*L/PFA 
1)2 (U * 02 (L ) ♦F»RML/PFA 
0(L) = Rl*(Cl(L)+D2d) ) 

IF(0(L )-OP(L) )62,62.61 
61 TERM (I»J)=TM(t) 

DP(L) = OIL) 

IF(DP(L)-70C•162,999,999 
6 2 CONTINU E 
GC TG 6G 

999 OCP(I,J,L) =0PJ L) 

1C1 CONTINUE 
3C5 CONTINUE 

OC 90 J*1,1SL 
_ DC. 90 J = Jjl5l 

or 90 l=T,Cmax 

OCP11,JfL) = OOP(J,I,L) 

TFRM(I,J) = TFRM(J,I) 

9C CONTINUE 

OC 91 1=1,ISL 

_DC -91 0=1,1 SL . .. 

DC 91 1=1,UMAX 
__ . JJ=(J»LMAX)-dMAX-L) 

P~ROB=(CCP( I »J»L)-20C. )/50C. 

_ I F.1 PRC £19 ’,9 A, 9 A 

93 NPRCB=C 

GC TO 95 _ 

9A IF (PRC8-1. )9A2»9A1,9A1 

SECTION I - CALCULATION OF MOVEMENT TIMES 


B-ll - 


?AJ NPRCB=1000 









GO TO <35 

942 NPRCB»PRQB»10CC. _ 

95 WRITE OUTPUT TAPE 6, 709 , l‘,*J, L« J j'.NPROB , OOP Cl ,~J , L 5 , TERM II , J ) 
WRITE OUTPUT TAPE 9 , 720 ,I,JJ,NPROR 


709 

FORMAT 11CH 

415,110,2F1C.3) 

72C 

FORMAT J10HB 

15,15,110) 

91 

CONTINUE 



_WRITE OUTPUT TA PE 6*710 

WRITE OUTPUT TAPE 9,710 
71C FCRHATIIOHR • ) 

N* 1 

DO 75C W-1»ISL 
IF(IGNORE)7101,71C 1,7102 

.. 7.1C 1 . continue__ 

READ INPUT TAPES»75i^PIM ),S(N),S(N+l),S(N+2) ,5 IN*31» SI N*4), 
1 SIN + 5 ) , S(N+6)* S C N♦7)» $ IN* C)» S I N*9) 

7102 CONTINUE 
N*N-»LMAX 
KP=P(M) 

_ . WR I TE OUTP UT T APE 6*75,2,M,KP 

WRITE OUTPUT TAPE 9,752,M»KP 


751 

FORMAT IF1C.C,1CF5 

.0) 

752 

FORMAT(15HS 

15,110 

75C 

CONTINUE 



WRITE OUTPUT TAPE 

6,753 


WRITE OUTPUT TAPE 

9,753 

' 753* 

FCRMATIllHS 

• ) 


K = LMAX * ISL 



DC 754 N=1,K 



KS*XFIXF(S(N)) 



WRITE OUTPUT TAPE 

6,755,N,KS 


WRITE OUTPUT TAPE 

9,755,N,KS 

"7*55 ~* 

FORMAT 115HD 

15,110 

754 

CONTINUE 



WRITE OUTPUT TAPE 

6,756 


WRITE OUTPUT TAPE 

9,756 

756 

FORMAT(11HD 

• ) 


T4*T44 1. 



II = II* ~ 



DO 759 I=1,ISL 
DC 759 J=1,ISL 
DC 759 L*1,LMAX 
OOP!I,J,L)*C. 

_ __ 159 . TERM (I,J)*0. 

.: .— IGNORE ■ 1 -- ----- 

GC TO 1 
714 CONTINUE 

WRITE OUTPUT TAPE 6,757 
WRITE OUTPUT TAPE 9,757 

_ _ 757_ FORMAT (17H _ENO_CF _PR0_8LEM I 

WRITE OUTPUT ~TAPE *6,758 ** 

WRITE OUTPUT TAPE 9,758 
756 FORMAT 11HC) 

ENO FILE 9 
REWINC 9 

i_SECT I CN I - CALCULATION OF MOVEMENT, TIWES 

__WRJ TE OUTPUT _TAP£ 4,U3 

713 FORMAT (18H PROBLEM COMPLETED ) 

__CALI EXIT .__... . . 
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